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We calculated the 4 He trimer and tetramer ground and excited states with the LM2M2 potential 
using our Gaussian expansion method (GEM) for ab initio variational calculations of few-body 
systems. The method has extensively been used for a variety of three-, four- and five-body systems 
in nuclear physics and exotic atomic/molecular physics. The trimer (tetramer) wave function is 
expanded in terms of symmetric three-(four-)body Gaussian basis functions, ranging from very 
compact to very diffuse, without assuming any pair correlation function. Calculated results of the 
trimer ground and excited states are in excellent agreement with the literature. Binding energies 
of the tetramer ground and excited states are obtained to be 558.98 mK and 127.33 mK (0.93 mK 
below the trimer ground state), respectively. We found that precisely the same shape of the short- 
range correlation (tij < 4A) in the dimer appear in the ground and excited states of trimer and 
tetramer. The overlap function between the trimer excited state and the dimer and that between 
the tetramer excited state and the trimer ground state are almost proportional to the dimer wave 
function in the asymptotic region (upto~ 1000 A). Also the pair correlation functions of trimer and 
tetramer excited states are almost proportional to the squared dimer wave function. We then come 
to propose a model which predicts the binding energy of the first excited state of 4 Hejv (N > 3) 
measured from the 4 Hejv-i ground state to be nearly 2 {n-i) ^ 2 usm S the dimer binding energy B^. 



I. INTRODUCTION 

In early 1970's, Efimov pointed out a possibility of 
having an infinite number of three-body bound states 
even when none exists in the separate two-body subsys- 
tems OH!]- This occurs when the two-body scattering 
length is much larger than the range of the two-body 
interaction. As a candidate of such three-body states, 
Efimov discussed about the famous Hoyle state [IJ (the 
second + state at 7.65MeV in the 12 C nucleus) taking a 
model of three a particles (clusters of three 4 He nuclei) as 
well as about the three- nucleon bound state ( 3 H nuclei). 
In nuclear systems, the Borromean states, weakly bound 
three-body states though having no bound two-body sub- 
systems, are familiar but not classified as Efimov states. 

In atomic systems, triatomic 4 He (trimer) have been 
expected to have bound states of Efimov type since the 
realistic 4 He- 4 He interactions give a large 4 He- 

4 He scattering length (~ 115A), much greater than the 
potential range (~ 10A), and a very small 4 He dimer 
binding energy (~ 1.3 mK). (Experimentally, Ref. [Toj 
evaluated a scattering length of 104j^ g A and a binding 
energy of l.l+g'J mK). 

As is mentioned in recent reviews about the 4 He 
trimer [ll|, E3 (further references therein), i) a lot of 
three-body calculations using the realistic pair poten- 
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tials have shown that the 4 He trimer possesses two bound 
states with binding energies of nearly 126.4 mK and 2.3 
mK, ii) it is already rather well established that, if the 
4 He trimer excited state exist, it should be Efimov na- 
ture, and iii) it is suggested that the 4 He trimer ground 
state may be considered as an Efimov state since the 
ground- and excited-state binding energies move along 
the same universal scaling curve under any small defor- 
mation of the two-body potential (for details, see, e.g., 
Sec. Ill of Ref. [13]). Experimentally, the 4 He trimer 
ground state has been observed in Ref. [l4| to have the 
4 He- 4 He bond length of ll^g A in agreement with the- 
oretical predictions, whereas a reliable experimental evi- 
dence for the 4 He trimer excited state is still missing. 

Only very recently, experimental evidences of Efi- 
mov trimer states have been reported in the work us- 
ing the ultracold gases of cesium atoms [HI [Hj], pota- 
sium ato ms fl7j . lithium-7 atoms [HI, EH, an d lithium-6 
atoms 24J , in which the two-body interaction between 
those alkali atoms was manipulated so as to tune the 
scattering length to values significantly greater than the 
potential range. These experiments have been access- 
ing the study of a wide variety of interesting physical 
systems in the atomic and nuclear fields. Recently, the 
study extends to the Efimov physics and universality of 
four-atomic systems (tetramers). 

Though the interactions between 4 He atoms can not 
be manipulated, the study of 4 He trimer using the re- 
alistic pair potentials has been providing fundamental 
information to the Efimov physics. Now it is one of the 
challenging subjects to precisely investigate the structure 
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of 4 He tetramer using the realistic 4 He- 4 He potential. 

So far there exist in the literature a large number of 
4 He trimer calculations [25M39I ] giving well converged re- 
sults with the realistic 4 He- 4 He interactions. However, 
calculations of the tetramer remain limited [25|, I28U3T] ; in 
those papers, the binding energy of the tetramer ground 
state agrees well with each other, while that of the loosely 
bound excited state differs significantly from one another. 

Thus the main purpose of the present paper is to per- 
form accurate calculations of the 4 He tetramer ground 
and excited states using a realistic 4 Hc- 4 Hc interaction, 
the LM2M2 potential [6j. We employ the Gaussian ex- 
pansion method (GEM) for ab initio variational calcu- 
lations of few-body systems (40lj43| . The method has 
been proposed and developed by the present authors 
and collaborators and applied to various types of three-, 
four- and five-body systems in nuclear physics and exotic 
atomic/molecular physics (cf. review papers (43l - l45| ). 

Advantage of using the GEM for the 4 He tetramer 
calculation in the presence of the strong short-range re- 
pulsive potential is as follows: Some 30000 symmetrized 
four-body Gaussian basis functions, ranging from very 
compact to very diffuse, are constructed on the full 18 
sets of Jacobi coordinates without assuming any pair cor- 
relation function. They forms a nearly complete set in 
the finite coordinate space concerned, so that one can de- 
scribe accurately both the short-range structure and the 
long-range asymptotic behavior (up to ~ 1000 A) of the 
four-body wave function, which makes it possible to find 
new facets of 4 He clusters. 

We thus find that precisely the same shape of the short- 
range correlation (r*jj < 4 A) in dimer appears in the 
ground and excited states of trimer and tetramer. This 
gives a foundation to an a priori assumption that a two- 
particle correlation function (such as the Jastrow's) so 
as to simulate the short-range part of the dimer wave 
function is incorporated in the trimer and tetramer wave 
functions from the beginning. 

By illustrating the asymptotic behavior of the 4 He 
trimer and tetramer, we discuss about an interesting 
relation between their excited-state wave functions and 
the dimer wave function. We then come to propose a 
'dimerlike-pair' model that predicts the binding energy 
of the first excited state of the iV-cluster system, 4 He at, 
measured from the ground state of 4 HeAr_i to be approx- 
imately 2 (N-i) ^ 2 usm S the dimer binding energy B 2 - 

We explicitly write the asymptotic form of the total 
wave function of 4 He trimer (tetramer). The asymptotic 
normalization coefficient (ANC) [H, EH, EHH , namely 
the amplitude of tail function of the dimer-atom (trimer- 
atom) relative motion in the present case, is a quantity to 
reflect the internal structure of trimer (tetramer) . There- 
fore, attention to the ANC might be useful when one 
intends to reproduce the non-universal variation of the 
He trimer (tetramer) states by means of parametrizing 
effective models beyond Efimov's universal theory. 

The paper is organized as follows: In Sec. II, we apply 
the GEM to the three-body calculation of the 4 He trimer 



ground and excited states showing that the calculated 
results agree excellently with the literature. In Sec. Ill, 
the four-body calculation of the 4 He tetramer ground and 
excited states is presented. Summary is given in Sec. IV. 

II. 4 He TRIMER 

The 4 He trimer bound states have extensively been 
studied in many theoretical work using realistic poten- 
tials. Monte-Carlo, hyperspherical, variational and Fad- 
deev techniques were used to calculate accurately the 
binding energies of the ground and excited states [25l — 
|39| (see also recent reviews [III [Hj]). Nevertheless, in 
this section, we explain our Gaussian expansion method 
(GEM) and present the calculated result for the 4 He 
trimer in order to demonstrate high accuracy of our cal- 
culation before we report our investigation of the He 
tetramer in the next section. 

A. Three-body wave function 

We take all the three sets of Jacobi coordinates (Fig.[l}, 
xi = r 2 — r 3 and yi = ri — i(r 2 + r 3 ) and cyclically for 
( X 2,Y2) and (x3,y3), being the position vector of ith 
particle. Hamiltonian of the system is expressed as 

b --&-& + £ v ™ (2 - i> 

where fi x = \m and fi y = |rn, m being mass of a 4 He 
atom. V(rij) is the two-body 4 He- 4 He potential as a 
function of the pair distance = r j — r s ; . 




FIG. 1: Three sets of the Jacobi coordinates for 4 He trimer. 

We calculate the three-body bound-state wave func- 
tion, ^3, which satisfies the Schrodinger equation 

(JJ- £0*3 = 0. (2.2) 

Since we consider the 4 He atom as a spinless boson, 
we expand the wave function of three identical spinless 
bosons in terms of L 2 -integrable, fully symmetric three- 
body basis functions: 

Q max 

VF3 = £ A a <P^ m \ (2.3) 

a=l 

$(sym) = $ Q ( Xl , yi ) + $ Q ( X2 , y2 ) + $ Q (x3,y 3 ). (2.4) 
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It is of importance that those basis functions {<I>Q yin - ) ; a = 
1, ...,a max }, which are nonorthogonal to each other, are 
constructed on the full three sets of Jacobi coordinates; 
this makes the function space of {$i sym) } quite wide. 

The eigenenergies E and amplitudes A a of the ground 
and excited states are determined by the Rayleigh-Ritz 
variational principle: 



H - E I # 3 ) = 0, 



(2.5) 



where a = l,...,a max . Eqs.(2.5) results in a generalized 
eigenvalue problem: 



£ [Hca- - EAf a ^,] A a , = 0. 

q' = 1 

The matrix elements are written as 

Ha, a' = ( $i Sym) | H I $^ m) ; 

M ataf = ($i sym) | l |*i? m) 



(2.6) 



(2.7) 
(2.8) 



The lowest-lying two 5*-wave eigenstates, ^^\v = 0, 1), 
will be identified as the trimer ground (v = 0) and excited 
(v = 1) states. 

We express each basis function $ Q (xi, y^) as a product 
of a function of Xj and that of y j : 



*a(x,:,y l ) = <t>n x i*(xi)il) n i (yt) Yi x (Ki)Yi (y*) 



where a specifies a set of quantum numbers 



JM 

(2.9) 



a = {n x l x ,n y l y , JM} 

commonly for the components i = 1,2,3. J is the total 
angular momentum and M is its z-component. In this 
paper, we consider the trimer bound states with J = 0. 
Then, the totally symmetric three-body wave function 
requires l x = l y = even. 

One of the most important issues of the present vari- 
ational calculation is what type of radial shape we use 
for <j)n x l x {xi) and ^Pn y i y {Vi)- The basis functions should 
be capable of precisely describing the strong short-range 
correlation (without assuming any correlation function a 
priori) and the long-range asymptotic behavior of very 
loosely bound states. 

The GEM recommends two types of functions which 
are tractable in few-body calculations and work accu- 
lately. One is the Gaussian function and the other, 
more powerful one, is the complex-range Gaussian func- 
tion [43[ . In the next subsection, we introduce the former 
that was successfully used in our previous study (Sec. 3.1 
of Ref. [43j]) of the 4 He trimer ground and excited state 
with the use of the HFDHE2 potential @ . The latter 
function is introduced in Sec. II. C. 



B. Gaussian basis functions 

The radial function (f> n ~l x (%) in (2.9) is taken to be a 
Gaussian multiplied by x lnx (similarly for ifin y l y (y) )'■ 



(2.10) 
(2.11) 



where normalization constants are omitted for simplicity. 

Setting of the ranges by stochastic or random choice 
does not seem suitable for describing the strong short- 
range correlation and the long-range asymptotic behavior 
of the wave function. Any intended choice of the ranges 
is necessary. The GEM recommends to set them in a 
geometric progression: 



x n 



- x x a x * 



= 1, 

= 1, 



(2.12) 



■•A 



with common ratios a x > 1 and a y > 1. 



(2.13) 
This 

greatly reduces the nonlinear parameters to be opti- 
mized. We designate a set of the geometric sequence 
by {?i" lax , Xx, instead of {n" lax , xx, a x } and sim- 

ilarly for {fi.™ ax , yx, y n max } ; which is more convenient to 
consider the spatial distribution of the basis set. Opti- 
mization of the nonlinear range parameters is in principle 
by trial and error procedure but much of experiences and 
systematics have been accumulated in many studies using 
the GEM. 

The basis functions {<f> n i} have the following proper- 
ties: i) They range from very compact to very diffuse, 
more densely in the inner region than in the outer one. 
While the basis functions with small ranges are respon- 
sible for describing the short-range structure of the sys- 
tem, the basis with longest-range parameters are for the 
asymptotic behavior, ii) They, being multiplied by nor- 
malization constants for I <b n i) = 1, have a relation 



On l 



>n+kl) = 



2a k 



1 



-,2k 



Z+3/2 



(2.14) 



which tells that the overlap with the fc-th neighbor is 
independent of n, decreasing gradually with increasing k. 

We then expect that the coupling among the whole 
basis functions take place smoothly and coherently so 
as to describe properly both the short-range structure 
and long-range asymptotic behavior simultaneously. We 
note that a single Gaussian decays quickly as x increases, 
but appropriate superposition of many Gaussians can 
decay even exponentially up to a sufficiently large dis- 



tance. A good example is shown in Fig. 3 of Ref. [43 1 
for the 4 He dimer wave function (with the HFDHE2 po- 
tential) that is accurate up to ~ 1000 A with the use of 
the nonlinear parameters {?i max = 60, x\ — 0.14 A and 
= 700 A } (the same-quality dimer wave function is 
seen in Fig. 2 below in Sec. II. D using the complex-range 
Gaussians with the LM2M2 potential). 

A lot of successful examples of the three- and four-body 
GEM calculations are shown in review papers [43 45] and 
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in papers of five-body calculations [49|,[5(|. The examples 
includes our previous calculation of the ground and ex- 
cited states of 4 He trimer using the HFDHE2 potential; 
the binding energies were in good agreement with those 
given by a Feddeev-equation calculation |33j . As for the 
trimer wave function, we showed, in Figs. 3, 13 and 14 in 
Ref . {43l] . that the strong short-range correlation (x < 4 
A) and asymptotic behavior (up to x ~ 1000 A) of the 
trimer ground and excited states are simultaneously well 
described. Also, the three-body basis functions (2.9)- 
(2.13) together with the LM2M2 potential were used re- 
cently by Naidon, Ueda and one of the present authors 
(E. H.) [5l[ to study the universality and the three-body 
parameter of 4 He trimers. 



C. Complex-range Gaussian basis functions 

Before we proceed to the calculation of the 4 He 
tetramer ground and excited states, we improve the 
Gaussian shape of the basis functions so as to have more 
sophisticated (but still tractable) radial dependence. We 
then test the new basis in the calculation of the trimer 
states below. 

In Ref. [43j |. we proposed to improve the Gaussian 
shape by introducing complex range instead of the real 
one: 



x i e -(i+iu)(x/x n y 



(2.15) 



where n — 1, ...,n™ ax and x n are given by (2.12). Using 
4>ni {%) , we construct two kinds of real basis functions: 

^° s \x) = x l e"^) 2 coau(x/x n ) 2 
4>^ n} (x) = x l e -^ x ^ 2 smu{x/x n f 



(2.16) 



(2.17) 



where we usually take ui = 1. The three-body basis func- 
tion $ a (xj,yi) in (2.9) is replaced by 
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(2.18) 



(2.19) 



where a specihes a set 

a= {'cos' or l sui\ij,n x l x ,nyly, JM}. 

The new basis {<^„"j^ (x)} apparently extend the func- 
tion space from the old ones (2.10) since they have the os- 
cillating components; see Sec. 2. 4 and Sec. 2. 5 of Ref. 43 1 
for some examples taking this advantage in calculations 
of highly vibrational excited states (with ~ 25 nodes) and 
scattering states. The sin-type basis (2.17) particularly 
work when the wave function is extremely suppressed at 
x ~ due to the strongly repulsive short-range potential. 

In the following calculations, we employ the new basis 
(2.16) and (2.17) for the a;-space instead of (2.10), but 
keep (2.11) for the j/-space. 



Note that, when calculating the matrix elements (2.7) 

and (2.8) using (x) and 0^ ln ' (x) , we explicitly take 
(2.15) and the right-most expression of (2.16) and (2.17) 
since the computation programming is almost the same 
as that for (2.10) though some of real variables are 
changed to complex ones. 

A great advantage of the real- and complex-range 
Gaussian basis functions is that the calculation of matrix 
elements (2.7) and (2.8) is easily performed. As for the 
overlap and kinetic-energy matrix elements of the trimer 
(tetramer), all the six- (nine-) dimensional integrals give 
analytical expression. In the case of the potential ma- 
trix, we have analytical expression except for the one- 
dimensional numerical integral having the final form 



x 2m e -\i 



V(x) x 2 dx . 



(2.20) 



We explained, in Ref. [43( , various techniques to perform 
the three- and four-body matrix-element calculations as 
easily, accurately and rapidly as possible. 

It is to be emphasized that the GEM few-body cal- 
culations need neither introduction of any a priori pair 
correlation function (such as the Jastrow function) nor 
separation of the coordinate space by x < r c and x > r c , 
r c being the radius of a strongly repulsive core potential. 
Proper short-range correlation and asymptotic behavior 
of the total wave function are automatically obtained by 
solving the Schrodinger equation (2.2) using the above 
basis functions for ab initio calculations. 



D. Pair interaction and He dimer 

To describe the interaction between the 4 He atoms, 
we employ one of the most widely used 4 He- 4 He inter- 
actions, the LM2M2 potential by Aziz and Slaman jg]. 
Use is made of — = 12.12 KA 2 as the input mass of 4 He 
atom. We can then precisely compare calculated results 
for the tetramer ground and excited states with those 
obtained by Lazauskas and Carbonell [25[ who made a 
Faddeev-Yakubovsky (FY) equation calculation taking 
the same potential and 4 He mass as above. Recently, the 
authors of Ref. [52j claim that a more precise value of 

= 12.11928 KA 2 should be employed. We shall ad- 
ditionally show the trimer and tetramer binding energies 
in the case of using this value. 

We calculated the 4 He dimer binding energy, say B2, 
and the wave function, \E , 2(= ^ / 2(a;)Yoo(x)), using the 
same prescription as described above. We expanded 
&2(x) with 100 basis functions of (2.16) and (2.17) as 

n max 

* 2 (x) = ]T [4 COB Vio S) (*) + ^ in ^ n) (*)] (2-21) 

71=1 

with a parameter set 

{n™ ax = 50, xi = 0.5 A, i„ r = 600.0 A, a = 1.0}. 

(2.22) 
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FIG. 2: Short-range structure (upper) and asymptotic be- 
havior (lower) of the radial wave function ^2(2:) of the 4 He 
dimer obtained by using the complex-range Gaussian basis 
functions (2.21) and (2.22). The open circles stands for the 
exact asymptotic form. The dotted line (upper) illustrates 
the LM2M2 potential in arbitrary unit. 



TABLE I: Mean values for 4 He trimer ground and excited 
states with the use of the LM2M2 potential and — = 12.12 
KA 2 . B ( 3 V) is the binding energy, ry stands for interparti- 
cle distance and j^g is the distance of a particle from the 
center-of-mass of the trimer. See text for the asymptotic nor- 
malization coefficient C3 (v = 0, 1). 



trimer 




ground state 






present 


Ref.[25j 


Ref.[26] 


Ref.[2JJ 


B^> (mK) 


126.40 


126.39 


126.4 


126.40 


(T) (mK) 


1660.4 


1658 


1660 




(V) (mK) 


-1786.8 


-1785 


-1787 




y/tfi) (A) 


10.96 


10.95 


10.96 




{r«) (A) 


9.616 


9.612 


9.610 




K 1 ) (A- 1 ) 


0.134 


0.135 






(r- 2 ) (A" 2 ) 


0.0228 


0.0230 






(A) 


6.326 




6.49 


6.32 


ci 0) (k-i) 


0.562 




0.567 




trimer 




excited state 






present 


Ref.[25] 


Ref.[26] 


Ref.[2JJ 


(mK) 


2.2706 


2.268 


2.265 


2.2707 


(T) (mK) 


122.15 


122.1 


121.9 




(V) (mK) 


-124.42 


-124.5 


-124.2 




VH> (A) 


104.5 


104.3 


101.9 




(m) (A) 


84.51 


83.53 


83.08 




<^> (A- 1 ) 


0.0265 


0.0267 






(r" 2 > (A" 2 ) 


0.00216 


0.00218 






(A) 


60.33 




58.8 


59.3 




0.179 




0.178 





We obtained B 2 = 1.30348 mK, y^F) = 70.93 A, and 
(x) = 52.00 A which arc the same as those obtained in 
the literature. Experimentally, (x) = 52 ± 4 A [l(| from 
which B2 = l.llg J m K was estimated. 

As shown in Fig. 2, both the strong short-range correla- 
tion (x < 4 A) and the asymptotic behavior of the dimer 
are well described. In the lower panel, i$ 2 W precisely 
reproduces the exact asymptotic shape 0.1498 exp(—K 2 a;) 
(A" 1 / 2 ) with k 2 = VmTh/h = 0.0104 A" 1 up to x-1200 
A which is large enough for our discussions. 

There are 30 basis functions whose Gaussian ranges 
x„ x < 4 A, which is sufficiently dense to describe the 
short-range structure of the wave function precisely. An 
interesting issue is whether the same shape of the short- 
range correlation in Fig. [5] appear also in the trimer and 
tetramer ground and excited states without assuming any 
two-body correlation function. 



E. Trimer bound states 

We calculated the wave functions of the trimer ground 
state, v t , 3°' ) , and the excited state, ^3 , and their bind- 
ing energies, B^ and Bg , respectively, as well as some 

mean values with the (v = 0, 1). Some of results are 
summarized in Table I together with those obtained in 
the literature. Our results excellently agree with those 
by Refs.[25U27j. The 4 He- 4 He bond length in the trimer 
ground state was measured as (r^) = A [T(ij |. which 
is well explained by the calculations, (ry) = 9.61 A. 

Those converged results were given by taking the 
symmetric three-body basis function {$Q Sym ^;a = 
l,...,a max } with a max = 4400, in which the shortest- 
range set is (xi=0.3A, yi =0.4 A) and the longest-range 
one is (x max = 150 A, y max = 600 A). All the nonlinear pa- 
rameters of the Gaussian basis set are listed in Table II. 
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TABLE II: All the nonlinear parameters of the Gaussian basis 
functions used for the 4 He-trimer states with J = (l x = ly). 



Those in column a) are commonly for 4>„i (x) and <f> * 1 " (x) 
and b) for tpn y l v (y). Total number of the basis, a max , is 4400. 



a) 


, (cos 


(x) 


J,(sin) 


(x) 




b) </> 


lyly 


(y) 




l x 




X\ 


Xn^}^ x 




ly 




in 




number 






[A] 


[A] 








[A] 


[A] 


of basis 





22 


0.3 


150.0 


1.0 





50 


0.4 


600.0 


2200 


2 


17 


0.6 


150.0 


1.0 


2 


40 


0.8 


400.0 


1360 


4 


14 


0.8 


130.0 


1.0 


4 


30 


1.0 


200.0 


840 



There are neither additional parameter nor assump- 
tions. The present calculation is so transparent that it 
is possible for the readers to repeat the calculation and 
check the results reported here. The parameters for the 
Gaussian ranges are in round numbers but further op- 
timization of them do not improve the binding energies 
(B^ 0) = 126.40 mK and B ( 3 1] = 2.2706 mK) as long as we 
calculate them with five significant figures (cf. another 
check in Sec. II. H about the accuracy of the calculation). 



Convergence of the binding energies and B^' with 
respect to increasing partial waves l x {= ly) is shown in 
Table III in comparison with the Faddeev calculation by 
Lazauskas and Carbonell [25j . The case l Xmxx — 4 is 
sufficient in the present work as long as the accuracy of 
five significant digits is required. 

The convergence of the present result is more rapid 
than that of the Faddeev solution (the same will be seen 
in the tetramer calculation in Sec. III). The reason is that 
both the interaction and the wave function are truncated 
in the angular-momentum space (/ max ) in the Faddeev 
calculations, but the full interaction is included in the 
present calculation (with no partial- wave decomposition) 
though the wave function is truncated (l m&K ). The differ- 
ence of the convergence in the two calculation methods 
was precisely discussed in the case of the three nucleon 
bound states ( 3 H and 3 He nuclei) in our GEM calcu- 
lation (4lTj43| and in a Faddeev calculation [53j : for an 
illustration of the difference, see Fig. 15 in Ref. (43|. In 
this context, it is worth pointing out that, in Table HI our 
result precisely agrees with another Faddeev calculation 
by Ref.[27] with no the partial wave decomposition. 



Use of the value ^- = 12.11928 KA 2 521 results in 

m 1 — IJ 

B^ 0) = 126.499 mK and B ( 3 1] = 2.27787 mK, while 
Ref. [HI gives 126.499 mK and 2.27844 mK, respectively. 
Calculation of the binding energy was also made pertur- 
batively with B 3 = 12 ^\ 9 2 28 ( T ) +{V), where (T) and 

(V ) are those obtained with ^- = 12.12 KA 2 ; this gives 
B^ = 126.498 mK and B^ = 2.27787 mK. The calcu- 
lations below in Sec.II.F-H take — = 12.12 KA 2 . 



TABLE III: Convergence of the 4 He trimer calculations 
with respect to the increasing maximum partial wave (£ ma x). 
The four columns present trimer ground (-B3 ') and excited 
(Sg 1 ') state energies in comparison with those obtained by 
the Faddeev-equation calculation of Ref. [25T| . 



trimer 

^max 


present 


Ref. [25] 


B< 0) (mK) B« (mK) 


s£ 0) (mK) B 


« (mK) 





121.00 2.2397 


89.01 


2.0093 


2 


126.39 2.2705 


120.67 


2.2298 


4 


126.40 2.2706 


125.48 


2.2622 


8 




126.34 


2.2677 


12 




126.39 


2.2680 


14 




126.39 


2.2680 



Short-range correlation and asymptotic 
behavior 



In order to see how the present method describes the 
short-range structure of trimer, we calculated the pair 
correlation function (pair distribution function or two- 



body density) (x) defined by 



p( v) (x 1 )Y 00 (5i 1 ) = 



(v). 



lyi j 



(2.23) 



where the symbol ( ) yi means the integration over yi 
only. This integration gives an analytical expression ow- 
ing to the use of the Gaussian basis functions; here, we 
explicitly rewrite as a function of (xi,yi) by trans- 
forming the other coordinates (x2,y2) and (x3,ya) into 
(xi,yx). (xi) is independent of i (— 1, 2, 3) and is ap- 
parently normalized as J P^ v \x)x 2 dx = 1. It presents 
the probability of finding two particles at an interparticle 
distance x. 

In Fig. [3l short-range structure of P^"\x)(v = 0, 1) is 
illustrated together with P 2 (x)(= \^ 2 (x)\ 2 ) for the 4 He 
dimer. The dashed line is for the trimer ground state 
(v = 0). The solid line for the excited state (v = 1) 
and the dotted line for the dimer have been multiplied 
by factors 14.5 and 6.0, respectively. It is of interest that 
precisely the same shape of the short-range correlation 
(x < 4 A) as seen in the dimer appears both in the trimer 
ground and excited states (the same will be seen in the 
tetramer ground and excited states in the next section). 
This gives a foundation to an a priori assumption that 
a two-particle correlation function (such as the Jastrow 
function) to simulate the short-range part of the dimer 
radial wave function ^2(x) is incorporated in the three- 
body wave function from the beginning. 

To investigate the trimer configuration in the asymp- 
totic region where one atom is far from the other two, 
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P^'(x) of the 



3: Short-range structure of the pair correlation function 
4 He trimer calculated by (2.23). The dashed 
line is for the trimer ground state (v = 0), the solid line for 
the excited state (v = 1) and the dotted line for the 4 He dimer 
(I^MI 2 )- The solid and dotted lines have been multiplied 
by factors 14.5 and 6.0, respectively, to be normalized at the 
peak. The same shape of the short-range correlation (x < 4 
A) appears in the three states. 



FIG. 4: Overlap function 0^\y), multiplied by y, between 
the 4 He trimer wave function (v = 0, 1) and the dimer one, 
which is denned by (2.24) . Open circles represent the fit of the 
asymptotic function (2.25) to 0^(y) using the asymptotic 
normalization coefficient Cg"K The solid line (v = 1) is found 
to be parallel to the dotted line for the dimer wave function 
(j/*2(y)). 



we calculate the overlap function 0^\y l ) [H EH, El-El] 
to describe the overlap between the trimer wave function 
^z\ v = 0: 1) an( l the dimer one ^(x): 

Oi" ) (yi)yoo(yi) = <*2(xi)|^ ) ) Xl . (2.24) 

In Fig. HJ we plot yO^\y) for the ground and excited 
states. They should asymptotically satisfy 

yO^iyf-^C^^i-K^y), (2.25) 

(v) 

where is the binding wave number given by 

4 V) = \/^ y (Bi v) - b 2 )/h (4 0) - o.ii7 a- 1 , 4 X) = 

0.0103A" 1 ). The amplitude C^ v) is called the asymptotic 
normalzation coefficient (ANC) [H, EJ [46M48j ] defining 
the amplitude of the tail of the radial overlap function. 
The asymptotic functions (2.25) with cf } = 0.562A"2 
and = 0.179A~5 (see the open circles) are precisely 
reproduced by the dashed line (v = 0) and the solid line 
(v = 1), respectively, up to y ~ 1000 A, which demon- 
strates the accuracy of our wave functions in the asymp- 

(v) 

totic region. The values of C3 agree with those given 
by Barletta and Kievsky [26j using a variational method 
with correlated hyperspherical harmonics functions (see 
Table I). 

The total three-body wave function ^^(v = 0, 1) is 
represented asymptotically as 

*W d v) ]T * 2 ( Xi )^^loo(yi). (2.26) 



The ANC, C3 , is a quantity to convey the interior struc- 
tural information of the trimer to the asymptotic behav- 
ior. It is known, in the nuclear peripheral reactions where 
only the asymptotic tails of the wave functions of react- 
ing particles contribute to the reaction process, the cross 
section is proportional to the squared ANC which can 
be measured in some specific systems [53, EH, HH ■ The 
idea of ANC might be available to the calculation of 4 He 
atoms reactions such as dimer + dimer — > trimer + atom. 
Also, attention to the ANC might be useful when one 
intends to reproduce the non-universal variations of the 
trimer states by parametrizing effective models. 

G. 'Dimerlike-pair' model in asymptotic region 

In Fig. El we find that the solid line (v — 1) is parallel 
to the dotted line (dimer); namely, k^\= 0.0103 A -1 ) 
is very close to «2(= 0.0104 A -1 ). This agreement is not 
accidental, but is understandable from a model, which 
we refer to as a 'dimerlike-pair' model, for the asymp- 
totic behavior of the trimer excited state (Fig. [S£i). The 
model tells that i) particle a, located far from b and c 
which are loosely bound (dimer) , is little affected by the 
interaction between b and c, ii) therefore, the pair a and 
b at a distance x is asymptotically dimerlike, iii) since 
x ~ y asymptotically, the amplitude of particle a along 
y is dimerlike, namely ~ k 2 . 

If this model is acceptable, we can predict that, in 
the asymptotic region, the pair correlation function of 
the trimer excited state, x 2 P^\x), should decay ex- 
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a) trimer excited state 
b 



dimerlike pair 
x 



b) tetramer excited state 

k dimerlike pair 



0^- 



OJ 



FIG. 5: 'Dimerlike-pair' model for the asymptotic behavior 
of the trimer and tetramer excited states (see text). 



ponentially with the same rate as that in the dimer 
(x 2 P2(x)). This is clearly seen in Fig. [HI the solid and 
dotted lines have almost the same exponentially-decaying 

rate of 2/4^ 2k 2 ). The same evidence is seen in Figs. 3 
and 14 of our previous calculation of the 4 He trimer using 
the HFDHE2 potential reported in Ref. [Hj]. 



Once we accept the dimerlike-pair model (k 



(i) 



to), 



we can estimate B 3 , the trimer excited-state binding 
energy, using B 2 . With the use of the definitions of the 
binding wave numbers: 



K 2 = y/2fl x B 2 /h, 



.(1) 



where /j, x = and fi y 
then predict 



2tLy(B™ - B 2 )/h, 



|m. Taking k 3 



(i) 



B 2 + ^B 2 = -B 2 = 2.281 mK, 

fiy 4 



(2.27) 
(2.28) 

k 2 , we can 
(2.29) 



which is close to 2.2706 mK by the present three-body 
calculation using the LM2M2 potential for which we have 



the ratio B 3 ' /B 2 = 1.74 (~ 7/4). 



In order to see a deviation of the ratio from 7/4 de- 
pending on the realistic potentials in the literature, we 
refer to B { 3 ] /B 2 = 1.59 (SAPT2 @), 1.65 jH, H3| 
(SAPT2007 (|), 1.74 [Hlg (TTY Q), 1.74 (LM2M2), 
2.01 [13] (HFDHE2). The dimerlike-pair model provides 
a reason why the ratio B 3 ^ /B 2 is located around 7/4 in 
a narrow region of 1.6-2.0. 

We note that this model should be considered under 
the condition that the 4 He atoms are interacting with 
a realistic pair potential and should not be discussed in 
any situation where a large deformation of the strength 
is posed to the potential (cf. a discussion in Sec. Ill of 
Ref. [l3[ on the Efimov states in 4 He trimer). 



He trimer 




500 

X (A) 



1000 



FIG. 6: Asymptotic behavior of the pair correlation (distri- 
bution) function P3 (sc), multiplied by x 2 , of the 4 He trimer. 
The dashed line is for the trimer ground state (v = 0). The 
solid line for the excited state (v = 1) and the dotted line for 
the dimer are found to have almost the same exponentially- 
decaying rate, 2k^ ~ 2k2- 



We try to apply the same model to the tetramer 
excited state (Fig. \5jp) and predict its binding energy 

B^K Asymptotically, particle a decays from the trimer 
(i). 



(b + c + d) as exp(— k\ z) with 



n^ = J2„ z (B^-B^)/n, 



?(!) 



(2.30) 



where fi z = |m is the reduce atom-trimer mass. Taking 



K%, we predict B^ as 



B 



(i) 



? (0) 



B 3 -I B 2 - B 3 



(0) 



\B 2 = 127.27 mK (2.31) 



when employing the calculated values of B 3 °^ and B 2 
with LM2M2. In Sec. Ill, we make a four-body calcula- 
tion of the tetramer with LM2M2 and check the above 



(i) 



prediction of B 



H. Generalized eigenvalue problem 

In this subsection, we discuss about a technical sub- 
ject on a numerical trouble which arises when solving 
the generalized eigenvalue problem (2.6). This is due to 
the fact that the overlap matrix Af becomes almost sin- 
gular when a very large number of nonorthogonal basis 

functions {$a ym ' 1 } employed. In this case, because of the 
non- negligible round-off error in double-precision compu- 
tation (~16 decimal digits), we may obtain no solution 
of (2.6) or a solution that includes some unphysically too 
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deep erroneous bound state. In order to overcome this 
trouble, we took the following two steps: 

Step i): we first diagonalize the overlap matrix Af: 



a' = l 



,n 



(AO 



(2.32) 



where a,N = l, a max . The eigenvalues z/jy are positive 
definite since Af ata > = Af a ',a- We then define a new, 
symmetrized orthonormal basis set: 



(sym) 
AT 



1 



($ (sy m) ! $ (sym) } 



3 AT, AT' 



(2.33) 
(2.34) 



where N, N' — l,...,a max . The generalized eigenvalue 
problem (2.6) are then equivalently converted into a stan- 
dard eigenvalue problem: 



TABLE IV: Stability of the calculated trimer binding ener- 
gies of the lowest-lying two states against decreasing number 
(Mnax) of symmetrized orthonormal three-body basis func- 
tions {3>jv ym ' ) } of Eq.(2.33). This assures accuracy of the con- 



clusion B { 3 0) = 126.40 mK and B^' = 2.2706 mK in Table I. 



(AN n - 



Bi 0) (mK) 



o(l) 
°3 



(mK) 



3250 






126.3999 


2.270606 


3240 


( 


-10) 


126.3998 


2.270605 


3200 


( 


-50) 


126.3995 


2.270602 


3150 


(- 


-100) 


126.3991 


2.270594 


2950 


(- 


-300) 


126.3975 


2.270533 


2750 


(- 


-500) 


126.3954 


2.270484 


2250 


(- 


1000) 


126.3657 


2.270163 



E 



N,N' 



E6 N:N ,] A N , = 0, 



N' 



where N = 1, ...,a max , and 



;( s y m ) 



(sym) 



(2.35) 



(2.36) 



Here, we arrange {® N } in the decreasing order of 



> vn > 



(2.37) 



When the nonorthogonality among the basis functions 
{$a y } is very large, some of vn become extremely 
small and therefore the large factor 1/sJvn may cause 
a serious cancellation in the summation in (2.33). Since 
the present calculation is performed by double-precision 
computation, such a large cancellation may generate a 
substantial round-off error in (2.33) and hence in the 
matrix elements (2.36). This may give rise to some erro- 
neous eigenstates in (2.35) that have unphysically huge 
binding energies. 

Step ii): We therefore omit such members of { < I ) ^ ym ''} 
that have too small Vfj. The binding energies of such 
unphysical states decreases quickly as the basis size is re- 
duced. Finally, we reach an appropriate size, say iV max , 
of the basis {$# y } f° r which those unphysical states 
have disappeared from the low-energy region, and ener- 
gies of the lowest-lying (deepest) states take physically 
reasonable values. It is to be emphasized that the bind- 
ing energies of so-obtained lowest-lying physical states 
are stable against further reduction of 7V max . 

Table IV explicitly demonstrates Step ii). We start 
with a max = 4400 basis functions {$« y } whose pa- 
rameters are given in Table II. When the size of the new 

basis { < I > ^ ym '} is reduced from a max to iV max = 3250 
according to (2.37), the solution of (2.35) has come to 
include no unphysical state and give the binging ener- 



gies B 



(0) 



126.3999 mK and B 



(i) 



the lowest two states. By checking the stability of the 
energy values against further decreasing A^ max , we verify 
the values of B ( ° ] = 126.40 mK and B^ = 2.2706 mK 
in Table I. 



III. 4 He TETRAMER 

Calculation of the 4 He tetramer using realistic poten- 
tials has been performed in Refs. p5l. [28l - l3lj |. Although 
binding energy of the ground state obtained in the pa- 
pers agrees well with each other (~ 558 mK), that of 
the loosely bound excited state differs significantly from 
each other; namely, the binding energy with respect to 
the trimer ground state (126.4 mK) is given as 1.1 mK 
by the Faddeev-Yakubovski (FY) equations method [25[ , 
6.6 mK by Monte Carlo methods combined with the adi- 
abatic hyperspherical approximation [30] and 52 mK re- 
cently by using a method_of the correlated potential har- 
monic basis functions 



31] 



2.270606 mK for 



Though the Faddeev result 
(1.1 mK) seems to the present authors the most accu- 
rate, the excited state was not solved as a bound-state 
problem in Ref. 25] but the result was extrapolated from 
the atom-trimer scattering phase shifts. 

Thus the purpose of this section is to perform, us- 
ing the same LM2M2 potential as in Ref. [25|], accurate 
bound-state calculation of the tetramer excited state, not 
only giving a precise binding energy but also describing 
the short-range correlation and the asymptotic behavior 
of the wave function properly. 

The GEM has extensively been employed in bound- 
state calculations of various four-body systems in nu- 
clear and hypernuclear physics (cf. review papers (43h - 
l^). Extension from three-body GEM calculations to 
four-body ones in the presence of strong short-range re- 
pulsion is a familiar subject in nuclear physics. For exam- 
ple, the study of three-nucleon bound states ( 3 H and 3 He 
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nuclei) in Ref. |4l[ was extended to that of four-nucleon 
ground state ( 4 He nucleus, J 71 " = + ) 15511 and the first 
excited, very diffuse state (J 71 ' = + ) [5a|. The study 
of the three-a-particle system ( 12 C nucleus) (43l [57l. [58j| 
was extended to that of the four-a-particle system ( 16 
nucleus) [HH with the strongly repulsive Pauli-blocking 
projection operator on the a-a motion. Therefore, ex- 
tension of the 4 He trimer calculation to the tetramer one 
is straightforward on account of those experiences. 



A. Method 

We take two types of Jacobi coordinate sets, K-type 
and H-type (Fig. [7]). Namely, for K-type, X! = r 2 — r 1 , 
yi = r 3 - i(ri+r 2 ) and zi = r 4 - ±(ri+r 2 +r 3 ) and cycli- 
cally for {xi,yi,Zi; i = 2,..., 12} by the symmetrization 
between the four particles. For H-type, X13 = r 2 — ri, 
yi3 = r 4 -r 3 , Z13 = i(r 3 +r 4 )- |(ri + r 2 ) and cyclically 
for {xi,yi,z;; i = 14, ...,18}. An explicit illustration of 
the totally 18 sets of the rearrangement Jacobi coordi- 
nates of four-body systems is seen in Fig. 18 of Ref. p3l |. 



with 




! 



'13 



o 

V13 



16 



63 



i=1 

K-type (/ =1-12) 



i=13 
H-type (/ =13-18) 



FIG. 7: Jacobi coordinates, K-type and H-type, for the 4 He 
tetramer. Symmetrization of the four particles generates the 
sets i = l, 12 (K-type) and i — 13, 18 (H-type). 

The total four-body wave function ^4 is to be obtained 
by solving the Schodinger equation 

(H-E)ty± = (3.1) 

with the Hamiltonian 

where fi x — \m, /j, y = |m and fi z = jin on the K-type 
coordinates, and /x x = [i y = and \i z — m on the 
H-type ones. 

$4 is expanded in terms of the symmetrized L 2 - 
integrable K-type and H-type four-body basis functions: 



$ ( r ;K) = f^O^y^), 

i=l 
18 

$ (sym;H) = V"^ 

CtH / j *BH 

i=13 



$W(x 4 , yi , Zl ), 



(3.4) 
(3.5) 



in which $(xi, y^, Zj) is a function of i-th set of Jacobi co- 



ordinates. It is of importance that $ 



(sym;K) 



and $ 



(sym;H) 



are constructed on the full 18 sets of Jacobi coordinates; 
this makes the function space of the basis quite wide. 

The eigenenergies E and amplitudes A^(A^) are 
determined by the Rayleigh-Ritz variational principle: 



( <^ m;K) \H-E\* 4 )=0, 
( $i s r ;H) \H-E\$ 4 )=0, 



(3.6) 
(3.7) 



where qk = 1, ...,aj£ ax and «h = F ...,ag ax . This set of 
equations results in a generalized eigenvalue problem: 

(3.8) 

mi The matrix ele- 



V V \n {c ' c ' ) - em {c ' c ' ] 1 a (c ' } = 

c'=K,H a c ,=l 



where c = K,H and a c = 1, ...,a 
ments are given by 



(3.9) 
(3.10) 



Up to here is the most general way of variational calcula- 
tions for bound states of identical spinless four particles. 
We describe the basis function (3>oh ) m the form 

^(XuYhZi) = <t>nTlS X i) <t>nvh(Vi) fnj z ( z i) 

>i[[Y la ^)Y lv (yi)] A Y lm (^)] jM , 

(* = 1,..., 12) (3.11) 

Cfx.J.^,) = 4>ZlS X i) ^n v l v {yi) fn^lA z i) 

x[[Y l ^Y ly {yi)] A YiM)\ JM , 

(i = 13,..., 18) (3.12) 

where ax specihes a set 

«k = {'cos' or 'sin', oj, n x l x , n y l y , n z l z , A, JM}, (3.13) 

which is commonly for the components i = 1, 12; and 
similarly for «h commonly for i = 13, 18. 

Since we consider the case of J = in this paper, 
the totally symmetric four-body wave function requires 
i) lx = even, l y + l z = even and A = l z for the K-type 



$ 4 = V" A^^^'^ + V" A^^^™' 1 ^, ( 3 - 3 ) b£LSiS ai1 ^ ^ lx = 6Ven ' ly = 6Ven aild ^ = lz = ° Ven fOT 



OK = l 



QH = 1 



the H-type basis. 
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In (3.11) and (3.12), the radial functions are assumed, 
as in Sec. II, to be 



COS Uj(x /%n x ) 



)2 , (3.14) 
(3.15) 
(3.16) 



with geometric sequences of the Gaussian ranges: 



Uriy 



- xia x - 



(n x = 1 
(n y = 1, n. 



.,< ax ), (3.17) 



in ax \ 

y > 



(3.18) 
(3.19) 



In (3.15), the 'cos(sin)'-type function is not adopted for 
^riyiyiv) of the H-type basis though y is the distance 
between two particles. This is because the 'cos(sin)'-type 
basis for the cc-coordinate are applied to all the pairs of 
H-type by the symmetrization of the four particles. 

In the tetramer calculation, the total number, a max — 
a™ + q™, of the symmetrized four-body basis func- 
tions (3.4) and (3.5) amounts to cv max ~ 30000, ranging 
from very compact to very diffuse, to obtain a well con- 
verged solution. Since the nonorthogonality among those 
basis functions is too large to solve directly the general- 
ized eigenvalue problem (3.8), we take the same two-step 
method as described in Sec. II. H in the trimer calcula- 
tion. We finally solve the same type of standard eigen- 
value problem as (2.35) using the symmetric orthonormal 
four-body basis functions, {$^ ym ' ) ; iV = 1, N max }, in 
which the basis with too small have been omitted. 



TABLE V: Convergence of 4 He tetramer calculations with 



respect to the increasing maximum partial wave (l n 



B 



and B4 1 ' are the binding energies of the ground and excited 
states, respectively. Column 'K+H' is the results with both 
the K-type and H-type basis functions in (3.3) and column 
'(K)' is that with the K-type basis only. The excited state 



is unbound (-B4 < B^' — 126.40 mK) in the cases denoted 



by the symbol ' — '. For comparison, B^ obtained by the 
FY-equations calculation [25j] is listed in the last column. The 
LM2M2 potential is employed. 



tetramer 



present 



(mK) 
K+H ( K ) 



Ref.[25] 
B^ (mK) B| 0) (mK) 



K+H ( K 



500.71 (185.96) - ( 

558.29 (508.62) 127.24 ( 
558.98 (532.56) 127.33 ( 



348.8 
505.9 
548.6 
556.0 
557.7 



TABLE VI: Nonlinear parameters of the four-body Gaus- 
sian basis functions, (3.11)-(3.19), used for the 4 He tetramer 
ground and excited states with J = (A = l z ) in the case 
of imax = 2 with 23504 basis functions (cf. Table V). In the 
leftmost column, K(H) stands for the K-type (H-type) basis. 
We take uj = 1.0 in (3.14). 



B. Binding energy 



i( cos ) A Bin ) 



In the calculation of the 4 He-tetramer ground state 



^4 and the excited state ^\ L> , the converged result was 
obtained by employing the symmetric four-body basis 
functions of a max = 29056 with l max — 4. Table V shows 
the convergence of the binding energies of the tetramer 
ground (B^) and excited (B^) states with respect to 
increasing / max . Column 'K+H' is the result with both 
the K-type and H-type basis functions in (3.3), and col- 
umn '(K)' is that with the K-type basis only. 

Contribution from the K-type basis is dominant, but 
that from the H-type is sizable. Without the latter the 



>(i) 



? (0) 



excited state does not become bound (B\ ' < B 
126.40 mK) even for l max = 4. Since both type bases are 
not orthogonal to each other, the role of H-type one can 
be substituted in principle by the K-type one if a very 
large Z max is employed. But, this is not practical; use of 
both types of bases is essentially important. 

As seen in Table V, convegence of the binding energies 
with increasing Z max is more rapid in our calculation than 
that in the FY-equations calculation 25[. The reason 



of this difference in the conversion is the same as that 
mentioned in the trimer calculation (Sec. II. E). 
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6 


0.3 


6.0 


2 


6 


0.8 


6.0 


2 


8 


1.0 


6.0 


H 


2 


6 


0.6 


6.0 





6 


0.8 


6.0 


2 


8 


1.0 


6.0 


H 


2 


6 


0.6 


6.0 


2 


6 


0.8 


6.0 


2 


8 


1.0 


6.0 



In Table VI we list the nonlinear parameters of the 
four-body Gaussian basis functions, (3.11)-(3.19), in the 
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TABLE VII: Stability of the calculated binding energies of 
the 4 He tetramer ground and excited states with respect to 
the number (iV max ) of the symmetrized orthonormal four- 
body basis functions corresponding to { < l?^ ym ' > } in Eq.(2.33) 
in the case of / ma x = 4 with the LM2M2 potential, 
concludes Bf ' = 558.98 mK and fif ' = 127.33 K. 



This 



N 

1 v max 


(AA> max ) 


(mK) 


(mK) 


24800 






558.983 


127.326 


24790 


(• 


-10) 


558.981 


127.326 


24780 


(• 


-20) 


558.980 


127.326 


24750 


(• 


-50) 


558.977 


127.326 


24700 


(- 


-100) 


558.975 


127.325 


24300 


(- 


-500) 


558.954 


127.323 


23800 


(- 


1000) 


558.924 


127.320 



case of l max = 2 with 23504 basis functions (cf. Table V) 
to avoid too long listing for Z max = 4. The range parame- 
ters are given in round numbers but further optimization 
of them do not improve the binding energies as long as 
we calculate them with five significant figures. 

When solving the generalized eigenvalue problem (3.8), 
we take the same two-step method as mentioned in 
Sec. II. H. Table VII shows stability of the binding en- 
erges of the ground and excited states {B^\v = 0,1) 
against the decreasing number N max of the symmetrized 
orthonormal four-body basis functions corresponding to 

{8<T°} in Eq.(2.33). 

Calculated binding energies and some of mean values 
of the tetramer ground and excited states are summer- 
ized in Table VIII(a) in comparison with those obtained 
by Lazauskas and Carbonell [25j with the FY-equations 
method in which the excited state was not obtained by 
a direct bound-state calculation but the binding energy 
(127.5 mK) was extrapolated from the atom-trimer scat- 
tering calculations. Our result of = 127.33 mK, 

which is very closed to in Ref.[25j], confirms the ex- 
istence of the very shallow bound excited state ( J = 0) of 
the 4 He tetramer. The tetramer excited state is located 
only by 0.93 mK below the trimer ground state (126.40 
mK). This is analogous to that the trimer excited state 
lies by 0.967 mK below the dimer; a reason was explained 
in Sec. II. G by taking the dimerlike-pair model. 



TABLE VIII: (a) Mean values of the 4 He tetramer ground 
and excited states calculated by the present work and the FY- 
equations method [25| with the use of the LM2M2 potential. 
B4 is the binding energy, r»j stands for interparticle distance 
and Tto is the distance of a particle from the center-of-mass 
of the tetramer. (b) The other literature work [2^ - [3l| on the 
binding energies (see text). The values originally given in 
units of cm -1 are transformed in units of mK in the paren- 
theses. 



tetramer 



ground state 



excited state 





present 


Ref.[25j 


present 


Ref. [25] 


B 4 (mK) 


558.98 


557.7 


127.33 


127.5 


(T) (mK) 


4282.2 


4107 


1639.2 




(V) (mK) 


-4841.2 


-4665 


-1766.5 




V<r&) (A) 


8.43 


8.40 


54.5 


34.4 


<r«) (A) 


7.70 




35.8 




(r-- 1 ) (A" 1 ) 


0.155 




0.0792 




(rlj 2 ) (A" 2 ) 


0.0285 




0.0117 




vV 2 G > (A) 


5.16 




33.3 






2.1 




0.10 




(b) 


tetramer 


Ref. [28] 


Ref. [29] 


Ref. [30J 


Ref. [31J 


Bf ] (cm" 1 ) 


0.388(1) 


0.3886(1) 


0.387(1) 


0.388 


(mK) 


(558) 


(559.1) 


(557) 


(558) 


B^ (cm" 1 ) 






0.0922 


0.124 


(mK) 






(133) 


(178) 



Table VIII(b) lists B { " ] and B\ x> obtained in other 
literature papers by the Monte Carlo methods |28-|3CI 
and by using the correlated potential harmonic basis [3 If . 

All Bf ] values agree well with the results by the present of the at0 m-trimer distance z as an extension from (2.24): 



h 

ment for the small difference of — as used in Sec. II. E 

m 

gives Bf ] = 559.23 mK and = 127.42 mK. The 

calculations below in Sec.III.C take ^ = 12.12 KA 2 . 

C. Short-range correlation and asymptotic 
behavior 

Definition of the pair correlation function (2.23) for the 
trimer is extended to the tetramer states ^^\v = 0, 1): 

Pi^x) Voo(Si) =<*i B) |*i o) > yi «, (3-20) 
where ( ) yi . Zl means the integration over yi and Zj.. 

(v) 

The overlap function between a tetramer state ^4 
and a trimer one y { 3 3) (v 3 = 0,1) is defined as a function 



and FY-equations calculations, but B^ by Refs. [301131] 
deviate significantly from our and FY-equations results. 



O ( ^' v \z 1 )Y 00 (z 1 ) 



- 1 1 it» 



/xi.yi ■ 



(3.21) 



Use of ■ 
mK and B 



12.11928 KA 2 521 results in B 



(0) 



559.22 



(i) 



127.42 mK. The same perturbative treat- 



All the integrals in (3.20) and (3.21) give the analytical 
expression owing to the use of Gaussian basis functions; 
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0.008- 




X (A) 50 100 

z (A) 



FIG. 8: Short-range structure of the pair correlation function 
P± v \x) of the 4 He tetramer calculated by (3.20). The dashed 
line stands for the tetramer ground (v = 0) state and the solid 
line for the excited (v = 1) state. For the sake of comparison, 
additionally shown are the dotted line for the trimer ground 
state and the dash-dotted line for the trimer excited state. 
The solid, dotted and dash-dotted lines have been multiplied 
by factors 2.76, 1.36 and 19.8, respectively; the same shape 
of the short-range correlation (x < 4 A ) appears in all the 
states (cf . Fig. [3] for dimer) . 

is to be transformed to a function of (xi,yi,zi), 
and *^ 3) to that of (xi,yi). 

In Fig. [SJ we illustrate the short-range structure of the 
pair correlation functions P 4 (x) of the tetramer ground 
(v = 0) and excited (v = 1) states together with those of 
the trimer states. It is to be emphasized that the same 
shape of short-range correlation (x < 4 A) appears in all 
the states (cf. Fig. 3 for the dimer) without introducing 
any pair correlation function. 

The pair correlation functions (x) of the tetramer 
(v = 0, 1) take very small values in the strongly-repulsive 
potential region (x < 1.5 A) in Fig. [8] relative ratio of the 
values to the peak value is ^10 -6 . This ratio is to be com- 
pared with ^10~ 2 in the case of the four-nucleon bound 
state ( 4 He nucleus) calculated with a realistic nucleon- 
nucleon interaction with a strong short-range repulsion 
(the ratio is seen in Fig. 1 of Ref. [55] for the calcu- 
lated pair correlation function of the 4 He nucleus). Wc 
understand that so strong is the repulsive core of the 
atom-atom potential. 

In Fig. [9j we plot the overlap function (z), mul- 

tiplied by z, between the tetramer states (v = 0, 1) and 
the trimer states («3 = 0, 1) in the region z < 100 A. The 
two lines of v = 0, 1(^3 = 0) are to be compared with the 
result in Fig. 4 of Ref.[25[ of the FY-equations calcula- 
tion; the latter result represents the K-type FY compo- 
nents as a function of atom-trimer distance z. The two 
kinds of the results are resemble to each other though 



FIG. 9: Overlap function 0^ a,v) {z) in (3.21), multiplied by 
z, between the trimer state (113 = 0, 1) and the tetramer state 
(v — 0, 1) as a function of the atom-trimer distance z. 

they do not stand for the same quantity. As for the ex- 
cited state, the FY component is derived approximately 
by modifying the K-type FY amplitude of the zero-energy 
scattering [251 ]; the resulting amplitude is slightly more 
enhanced in the inner region than our overlap function 
of v = l(i>3 = 0). This is reflected in the r.m.s distance 

J (rfj) in Table VIII(a). In the plot of the overlap func- 
tions between the trimer excited state (v^ = 1) and the 
tetramer states, it is reasonably seen that O4 3 ~~ 1,v (z) is 
much smaller than O^ 3 ~ 0,v \z) (v = 0, 1) and decreases 
more rapidly along z. 

In Fig. EE0J zO[ V3=0 - v) (z) (v = 0, 1) are illustrated in 
the asymptotic region. They should satisfy 

zO^=°' v \z) ^ Cf ) exp(~4 w) z) (3.22) 

with k£° = ^2fi z (B ( 4 ] -B^)/h (4° } =0.231 A" 1 and 

k { P = 0.0107 A" 1 ). The dashed line (v = 0) and solid 
line (v — 1) reproduces the asymptotic functions (3.22) 

with the ANC Cf ) = 2.1 A-i and cf' = 0.10 A~ a (see 
the open circles), respectively, up to y ~ 1000 A. 

Asymptotically the tetramer is dissociated into the 
trimer ground state and a distant atom in the symmetric 
way between the four atoms (with negligible amount of 
the trimer excited state). Therefore, the tetramer wave 
function ^^\v = 0, 1) is represented asymptotically as 

*C»> _^ c*0 J2 *£> e -^l Yoo (z n ), (3.23) 

n=l Zn 

where the summation over n symmetrizes the four atoms; 
namely, the nth atom is isolated at z n from the trimer 
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FIG. 10: Asymptotic behavior of the overlap function 
q(v 3 ~o,v) multiplied by z, between the trimer ground state 
(«3 = 0) and the tetramer states (v = 0, 1). Open circles rep- 
resent the fit of the asymptotic function (3.23) to O^" 3 ~ ' v \z) 
using the asymptotic normalization coefficient 



FIG. 11: Asymptotic behavior of the pair correlation func- 
tion P^ v \x), multiplied by x 2 , of the 4 He tetramer calculated 
by (3.20). The dashed line stands for the tetramer ground 
(v — 0) state and the solid line for the excited (v = 1) state. 
The dotted line for the trimer excited state and the dotted 
line for the dimer are added for the sake of comparison. 



^3 n in which the nth atom is absent and the other three 
atoms are symmetrized. 

As for the tetramer's ANC, we can make the same 
comments as those for the trimer's ANC below Eq.(2.26). 



of the reduced mass of the dimerlike pair {\rn) to that 
of the 4 He- 4 HeAr_i system (^-m). But, as discussed 
in Sec. II. G, the factor might depend on the realistic po- 
tentials with a small deviation (roughly ±0.3). 



D. 'Dimerlike-pair' model in asymptotic region 



IV. SUMMARY 



In Fig. [TQl we note that the exponentially-decaying 
slope of the solid line («4 ) is very close to that in the 
trimer excited state («3 ) and that in the dimer («2)- 
This gives a support to the dimerlike-pair model for the 
tetramer excited state in the asymptotic region (Fig.[5]D). 

We can therefore predict that, in the asymptotic re- 
gion, the pair correlation function x 2 P± (x) should be 

proportional to exp(— 2 ). This is seen in Fig. 1111 
though the solid line (v = 1) is not so excellently straight 
as in the trimer excited state (dot-dashed line) due to the 
complexity of the four-body calculation. 

As mentioned in Sec. II. G, the dimerlike-pair model 
predicted ~ 127.27 mK. In the present calculation, 

we obtained = 127.33 mK. The prediction is very 
good as long as the LM2M2 potential is employed. 

Generally, in the case of 4 HeAr, the excited-state bind- 
ing energy calculated using realistic potentials may be 
expected as 



B 



(i) 

N 



B 



(0) 

N-l 



N 



2(N - 1) 



Bo 



(3.24) 



We have calculated the ground and excited states of 
4 He trimer and tetramer using the LM2M2 potential 
which has a strong short-range repulsive potential and is 
one of the most widely used 4 He- 4 He interactions. We 
employed the Gaussian expansion method (GEM) for 
ah initio variational calculations of few-body systems. 
The symmetrized three-(four-)body Gaussian basis func- 
tions, ranging from very compact to very diffuse, are con- 
structed on the full sets of possible Jacobi coordinates. 
Therefore, the basis set, spanning a wide function space, 
is suitable for describing both the short-range correlation 
(without assuming any pair correlation function) and the 
long-range asymptotic behavior of the trimer (tetramer) 
wave function as well as suitable for obtaining accurate 
binding energies. The main conclusions are summarized 
as follows: 

i) Calculated binding energies of the trimer ground and 
excited states, B^ and B^\ respectively, agree excel- 
lently with the literature (Table I); we have B^ 0} = 126.40 



mK and B 



(i) 



2.2706 mK. 



ii) As for the binding energies of the tetramer ground 
and excited states, we obtained Bf ] = 558.98 mK and 



where the factor multiplied to Bi comes from the ratio B 



(i) 



127.33 mK (situated only 0.93 mK below the 
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atom-trimer threshold) . The former is in good agreement 
with the literature calculations, while the latter supports 
the result of 127.5 mK by Ref. [25| differently from the 
other literature results (Table VIII). 

iii) We found that the strong short-range correlation 
{jij < 4 A) seen in the dimer appears also in the ground 
and excited states of the trimer and tetramer precisely in 
the same shape (Figs. 3 and 8). This gives a foundation to 
an a priori assumption that a pair correlation function to 
simulate the short-range part of the dimer wave function 
is incorporated in the three-(four-)body wave function 
from the beginning. 

iv) Illustrating the overlap function between the trimer 
excited state and the dimer (Og" -1 ' 1 ) and that between 
the tetramer excited state and the trimer ground state 

we found that those overlap functions are 
almost proportional to the dimer wave function in the 
asymptotic region up to ~ 1000 A (Figs. 4 and 10). Also 
it was found that the pair correlation functions of trimer 
and tetramer excited states (Pjj° and P^" -1 ^, respec- 
tively) are almost proportional to the squared dimer wave 
function in the asymptotic region (Figs. 6 and 11). We 
then came to propose a 'dimerlike pair' model (Fig. 5) 



that predicts the excited-state binding energy of 4 He at 
(N > 3) using Eq. (3.24). It will be of interest to ex- 
amine this model in the case of N > 5. A five-body 
calculation of the pentamer, 4 He5, is in progress. 

v) We calculated the asymptotic normalization coef- 
ficient (ANC) of the tail function of the dimer-atom 
(trimer-atom) relative motion in trimer (tetramer). This 
result may be available in peripheral reactions (insen- 
sitive to the interior of the system) including the He 
trimers and tetramers, in which the reaction cross sec- 
tion will be proportional to the squared ANC. The ANC 
is a quantity to convey the interior structural information 
to the asymptotic behavior. Therefore, attention to this 
quantity might be helpful when one tries to reproduce 
the non-universal variation of the 4 He trimer (tetramer) 
states by means of parametrizing effective models beyond 
Efimov's universal theory. 
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